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ABSTRACT 

We present spectroscopy of the eclipsing recurrent nova U Sco. The radial velocity 
semi-amplitude of the primary star was found to be Kw = 93 ± 10 kms~^ from the 
motion of the wings of the Hell A4686A emission line. By detecting weak absorption 
features from the secondary star, we find its radial velocity semi-amplitude to be i^/j = 
170 ± 10 kms~^. From these parameters, we obtain a mass of Mi = 1.55 ± Q.2AMq for 
the white dwarf primary star and a mass of M2 = 0.88 ± O.ITM© for the secondary 
star. The radius of the secondary is calculated to be R2 = 2.1 ± O.2i?0, confirming 
that it is evolved. The inclination of the system is calculated to be « = 82.7° ± 2.9°, 
consistent with the deep eclipse seen in the lightcurves. The helium emission lines are 
double-peaked, with the blue-shifted regions of the disc being eclipsed prior to the 
red-shifted regions, clearly indicating the presence of an accretion disc. The high mass 
of the white dwarf is consistent with the thermonuclear runaway model of recurrent 
nova outbursts, and confirms that U Sco is the best Type la supernova progenitor 
currently known. We predict that U Sco is likely to explode within ^ 700 000 years. 

Key words: accretion, accretion discs - binaries: eclipsing - binaries: spectroscopic 
- stars: individual: U Sco - novae, cataclysmic variables. 
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1 INTRODUCTION 



U Sco Be longs to a small class of cataclygmie variabl e s (CVs) 
known as recurrent novae (RNe) , which show repeated nova 
outbursts on timescales of decades. Of these, U Sco has the 
shortest recurrence period (~ 8y), with recorded outbursts 
in 1863, 1906, 1936, 1979, 1987 and 1999. It is likely that 
others have been missed due to its proximity to the ecliptic 
(4°). 

The RN class is further divided into three groups, de- 
pending on the nature of the secondary star (red giant, sub- 
giant or dw arf) , each of w hich have different eruption mech- 
anisms (see Warner 1995 & Webbink et al. 1987 for reviews) . 
U Sco is the on ly object of th e three in the sub-giant class 
which eclipses ( Schaefer 199C), and als o shows evidence of 
secondary absorption lines (Hanes 1985), making it the best 



away (TNR) model of classica l novae outbursts (e.g. Star 
rfield, Sparks & Truran 1985). In this model, the primary 



candidate for the determination of system parameters. 

The outburst mechanism for the U Sco sub-class of RNe 
is believed to be a modified version of the thermonuclear run- 
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star builds up a layer of material, accreted from the sec- 
ondary star, on its surface. The temperature and density 
at the base of the layer become sufficiently high for nu- 
clear reactions to begin. After ignition, the temperature rises 
rapidly and the reaction rates run away, until the radia- 
tion pressure becomes high enough to eject most of the ac- 
creted material. This pr ocess occurs when the accreted layer 
reaches a critical mass (Truran & Livio 1986), the value of 



which is a strongly decreasing function of the white dwarf 
mass. Therefore, to reduce the time interval between erup- 
tions to those observed in the RNe, either the mass accre- 
tion rate (M) or the white dwarf mass must be increased. 
For RNe, the white dwarf mass is the more important con- 
straint because there is an upper limit to M above which the 
degeneracy on the surface of the white dwarf weakens, and 
no powerful eruption can occur. For this model to account 
for the short times between outbursts as seen in RNe, the 
white dwarf must have a mass close to the Chandrasekhar 
limit. A tight constraint on the white dwarf mass in U Sco 
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is hence theoretically important, as it will provide a direct 
observational test of the TNR model for this group of RNe. 
The mass of the white dwarf in U Sco is also important 
in terms of binary evolut ion and the role of RN e as Type 
la supernova progenitors. Livio & Truran (1994) suggested, 
on the basis of the abundance determinations of the ejecta 
of RNe, that the mass of the envelope ejected during the 
outburst might be smaller than the amount of accreted ma- 
terial. In the case of high mass white dwarf systems, the 
primary could be pushed over the Chandrasekhar limit to 
produce a supernova explosion. A tight constraint on the 
white dwarf mass and the mass accretion rate allow s the 
time to supernova to be calculat ed, which according to ptar 



rfield, Sparks fc Shaviv (1988) 
years 



could be as little as 50 000 



There have been two previous at tempts to measure the 



mass of the white dwarf in U Sco. Johnston & Kulkarni 
(1992)Mattempted to use low resolution spectra to deter- 
mine the radial velocity semi-amplitude of both compo- 
nents by looking at the emission and absorption line shifts. 
They calculated a white dwarf mass of O.23-O.6OM0, a re- 
sult disputed by Duerbeck et al. (1993), who obtained 10 
spectra with poor phase coverage, but were able to con- 
clude a high primary mas s, although their errors were la rge 



(Ml = 1.16 ± O.69M0). ISchaefer fc Ringwald (1995) 



phased the radi al velocity measurements of Johnston & 
Kulkarni (1992) and Duerbeck et al. (1993) using a new 
orbital period. They concluded that due to inconsistent 7 
velocities, a phase shift in the emission lines and a large 
scatter in the radial velocity curves, no reliable white dwarf 
mass could be derived from the data. 

In addition to this, Hachisu et al. (2000a) sucessfuUy 
modelled the quiescent light curve of U Sco with a white 
dwarf mass of Mi ~ 1.37M0 , a secondary mass of M2 =0.8- 
2.OM0 and an orbital inclination of ~ 80°. Indirect support 
that the TNR model is the outburst mechanism in U Sco has 
also been found by the observation of luminous supersoft X- 
ray emission after the 1999 outburst (Kahabka et al. 1999). 
They also constrain the white dwarf to be massive from the 
temperature of the optically thick supersoft component. 

In this paper, we present new measurements of the ra- 
dial velocity semi-amplitude of the primary and secondary 
stars in U Sco and hence a new determination of the mass 
of the white dwarf in this system. 



2 OBSERVATIONS 

On the nights of 1999 April 15-19 we obtained 51 spectra 
of the recurrent nova U Sco with the 3.9m Anglo- Australian 
Telescope (AAT) in Siding Springs, Australia. We covered 
one complete orbit of U Sco during this time, including an 
eclipse - a full journal of observations is given in Table hi 
The exposures were all 1500 s with about 50 s dead-time for 
archiving of data. The setup comprised of the RGO Spec- 
trometer + 250mm camera, the TEK 1024 x 1024 CCD chip 
and the 1200V grating, which gave a wavelength coverage 
of approximately 4620A-5435A at 1.6A (95 kms~^) reso- 
lution. We also took spectra of 20 class IV and V spectral 
type templates ranging from FO - K2, and the flux standard 
LTT9239 (Hamuy et al. 1992). The 1.5 arcsec slit was orien- 
tated to cover a comparison star ~ 1 arcmin west of U Sco 
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Figure 1. R-band lightcur ve of the 1999 outburst of U Sco. The 
data have been taken from Matsunioto, Kato & Hachisu (2001) 



The open circle represents a measurement with an unrecorded 



Table 1. Journal of observations for U Sco. Orbital phase is cal- 
culated using the new ephemeris presented in this paper (equation 



UT Date 


UT 


UT 


No. of 


Phase 


Phase 




start 


end 


spectra 


start 


end 


1999 April 15 


14:57 


19:09 


8 


0.31 


0.44 


1999 April 16 


14:14 


19:21 


11 


0.09 


0.26 


1999 April 17 


14:21 


19:11 


10 


0.92 


1.07 


1999 April 18 


14:20 


19:16 


11 


0.72 


0.88 


1999 April 19 


14:23 


19:30 


11 


0.54 


0.70 



in order to correct for slit losses. Arc spectra were taken be- 
tween every U Sco exposure to calibrate instrumental flex- 
ure. The nights were all photometric, and the seeing varied 
between approximately 1.0 arcsec and 1.8 arcsec throughout 
the five nights. 

U Sco underwent its sixth recorded outburst on 1999 
February 25, reaching a maximum brightness of V = 7.6 
(Schmeer et al. 1999). The observations recorded in this pa- 
per were taken approximately 49 days after outburst maxi- 
mum, as can be seen from the lightcurve in Fig. hi, by which 
time the magnitude had decreased to approximately one 
magnitude above its quiescent level. Unfortunately, the sin- 
gle data point in Fig. nl corresponding to the dates of our 
observations has no recorded uncertainty. It is therefore im- 
possible to conclude whether we actually observed U Sco 
during an anomalously high state. What is certain, however, 
is that the additional shot noise in the continuum made it 
more difficult than we expected to detect absorption lines 
from the secondary star. The total fiux in each spectrum is 
plotted in Fig. as a function of Heliocentric Julian Date 
(HJD) following equation |l| (see Section ^j). The eclipse 
occurs on the third night of observation. There is also some 
evidence for a slight downward trend in brightness over the 
five nights. 
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Table 2. Times of mid-eclipse for U Sco according to Schaefer ( 
Ringwald (1995; SR95) and this paper. 



Cycle 


HJD at mid-eclipse 


O-C 


Reference 


(E) 


(2,440,000-1-) 


(sees) 







7717.5968 ±0.0139 


-1530.6 


SR95 


1 


7718.8412 ± 0.0069 


-334.1 


SR95 


4 


7722.5318 ±0.0139 


-425.4 


SR95 


5 


7723.7783 ±0.0139 


952.5 


SR95 


8 


7727.4699 ± 0.0139 


947.6 


SR95 


1168 


9154.9111 ±0.0104 


1003.5 


SR95 


1496 


9558.5240 ±0.0069 


293.1 


SR95 


1497 


9559.7497 ±0.0069 


-126.1 


SR95 


1500 


9563.4432 ± 0.0139 


33.2 


SR95 


2900 


11286.2143 ±0.0050 


-141.3 


This Paper 



HJD (2451280+) 



Figure 2. Light curve of U Sco during our observations. The 
solid curve is a parabolic fit to the eclipse minimum. 



3 DATA REDUCTION 

We debiased the data and corrected for the pixel-to-pixel 
variations with a tungsten lamp flat-field. The sky was sub- 
tracted by fitting second-order polynomials in the spatial di- 
rection to the sky regions on either side o f the object spectra. 
The data were then optimally extracted ( Home 1986| ) to give 
raw spectra of U Sco and the comparison star. Arc spectra 
were then extracted from the same locations on the detector 
as the targets. The wavelength scale for each spectrum was 
interpolated from the wavelength scales of two neighbouring 
arc spectra. The root-mean-square error in the fourth-order 
polynomial fits to the arc lines was ~ 0.025A. 

The next stage of the reduction process was to correct 
for the instrumental response and slit losses in order to ob- 
tain absolute fiuxes. A third-order spline fit to the contin- 
uum of the spectrophotometric standard star was used to 
remove the large-scale variations of instrumental response 
with wavelength. The slit-loss correction was performed by 
dividing the U Sco spectra by spline fits to the comparison 
star spectra, and then multipling the resulting spectra by a 
spline fit to a wide-slit comparison star spectrum. 



4 RESULTS 

4.1 Ephemeris 

The times of mid-eclipse were determined by fitting a 
parabola to the eclipse minimum (the solid curve in Fig, f/ 



A least-squares fit to t he 9 eclipse timings found in Schae- 



fer fc Ringwald (1995) and our eclipse yields the following 



ephemeris: 



T, 



mid — eclipse 



HJD 2 447 717.6145 
± 0.0044 



+ 1.230 5522 B 
± 0.000 0024 



(1) 



The differences between the observed and calculated 
times of mid-eclipse are given in Table y. We fin d no evi- 
dence for a non-zero value for P, in agreement with Schaefer 
fc Ringwald (1995)[ 



4.2 Average Spectrum 

The average spectrum of U Sco is displayed in Fig. |3| and in 
Table y we list fiuxes, equivalent widths and velocity widths 
of the most prominent lines measured from the average spec- 
trum. 

The spectrum is dominated by single-peaked emission 
lines between AA4870-5080A originating from the recent 
outburst. The identification of these nebular lines is com- 
plicated by their large widths, but the blend is likely to 
have contributions from NH, NHI, NV, OIII and OV, in- 
dicating a very high degree of excitation. Both Sekiguchi 
et al. (1988) and Munari et al. (1999) attribute the complex 
around A5015A to a blend of He I and Nil in the 1987 and 
1999 outbursts, respectively. Barlow et al. (1981) suggest 
in the 1979 outburst the presence of H/3 at A4861A, and 
strong OHI emission at AA4959/5007A. This is supported 
by comparison with nebulae around classical novae such 



as Nova Cyg ni 1992 (Moro-Martin, Garnavich & Noriega- 
Crespo 2001). There is also evidence of CIV emission at 
A4658A, in agreement with the 1987 outburst (Sekiguchi 
et al. 1988). The other strong features in the blend are 
likely to be due to NHI and NV, in agreement with Bar- 
low et al. (1981), and optical spectra of Nova GQ Muscae 
1983 (Pequignot et al. 1993). It is difficult to directly com- 
pare our spectra to the spectra of the 1979 outburst (Barlow 
et al. 1981) and the 1987 outburst (Sekiguchi et al. 1988) 
because the latter two datasets were taken much closer to 
the eruption. 

U Sco also shows strong, broad, double-peaked He II 
lines at A4686A and A5412A, typical of emissio n from an 
accretion di sc in a high inclination system 
Marsh 1986]). 



(e.g. 



Home & 



The secondary star features cannot clearly be seen in 
Fig. H, as the spectra have been averaged without being cor- 
rected for orbital motion. Consequently, any faint absorption 
lines will have been smeared out. 

However, when the orbital motion of the secondary 
star is removed, the MgI6 triplet (AA5167,5173,5184A) can 
clearly be identified (see Section 4/7 and Figure |l0| ). 
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Figure 3. The average spectrum of U Sco, uncorrected for orbital motion. 



5400 



Table 3. Fluxes and widths of prominent lines in U Sco, measured from the average spectrum. The wavelength range of the blend is 
AA4870-5080A. 



Line 



Flux 



EW 

A 



FWHM 
kms~^ 



FWZI 
kms~^ 



Hell A4686A 
Hell A5412A 
Blend 



1.345 ±0.007 
0.177 ±0.005 
11.95 ±0.02 



6.98 ±0.04 
1.18 ±0.03 
82.9 ±0.1 



700 ± 50 
700 ± 100 



1600 ± 100 
1500 ± 300 



4.3 Light Curves 

The contirmum light curve for U Sco was computed by sum- 
ming the flux in the wavelength range AA5080-5390A, which 
is devoid of emission lines. A polynomial fit to the contin- 
uum was then subtracted from the spectra and the emission 
line light curves were computed by summing the residual 
flux after continuum subtraction. 

The resulting light curves are plotted in Fig. M as a 
function of phase, following our new ephemeris. The contin- 
uum shows a symmetrical eclipse, with flickering through- 
out. There appears to be a dip in magnitude around phase 
0.5, perhaps indicative of a secondary eclipse, but there are 
too few data points to be conclusive. In addition, as a con- 
sequence of U Sco's long orbital period (1.23 days), the 
lightcurves represent data from five nights. Therefore, not 
all of the variability can be necessarily attributed to orbital 
effects. 

The eclipses of the Hell lines have a different shape from 



the continuum. Both lightcurves show a long egress, and 
a pre-eclipse increase in brightness. There is no evidence 
for significant variability in the lightcurve of the nebular 
lines (AA4870-5080A), supporting the idea that the lines 
are produced in a region external to the orbits of the stellar 
components. 



4.4 Trailed spectrum 

We subtracted polynomial fits to the continua from the spec- 
tra and then rebinned the spectra onto a constant velocity- 
interval scale centred on the rest wavelength of the lines. The 
rest wavelength of the nebular lines was taken as A4975A. 
The data were then phase binned into 40 bins, 36 of which 
were filled. Fig. B shows the trailed spectra of the Hell 
A4686A, Hell A5412A and nebular lines in U Sco. 

The Hell lines show two peaks, which vary sinusoidally 
with phase. This, and the rotational disturbance (where the 



© 0000 RAS, MNRAS 000, 000-000 



1.5 



0.5 



„ 2 

"b 

> 1-5 

W 
U 

o; 

E 0.5 


0.4 

\ 

y 0.3 

so 
u 

OJ 

^ 0.2 

o 

^ 0.1 



1.5 



0.5 



The recurrent nova U Sco 



C$»o<**> 



% J» 



Continuum 
XA50a0-5390A 






V * 



♦> 



Hell X46a6S 









I, t 



I** *^ tint 



He 



^f ^ 



tf 



* 



1* ,t A 



•.V.V'*Vv-' --^-V-> 



V rc.^'W"'*'- 



Nebular Lines 
XA4870-5080A 



-0.5 



0.5 
Orbital phase 



1.5 



Figure 4. Continuum and emission line lightcurvcs of U Sco. The open circles represent points where the real data (the closed circles) 
have been folded over. 



blue-shifted peak is eclipsed before the red-shifted peak) is 
evidence for origin in a high inclination accretion disc. There 
is also some evidence for an emission component moving 
from blue to red between phases 0.6 - 0.9. 



4.5 Doppler Tomography 

Doppler tomography is an indirect imaging technique which 
can be used to determine the velocity-space distribution of 
line emission in cataclysmic variables. For a detailed review 



of Doppler tomography, see Marsh (2000) 



Fig. M shows the Doppler map of the Hell A4686A line in 
U Sco, computed from the trailed spectra of Fig. fel but with 
the eclipse spectra removed. The data for the Hell A5412A 
line were too noisy to produce a Doppler map. The three 
crosses on the Doppler map represent the centre of mass 
of the secondary star (upper cross), the centre of mass of 
the system (middle cross) and the centre of mass of the 
white dwarf (lower cross). These crosses, the Roche lobe of 
the secondary star, and the predicted trajectory of the gas 
stream have been plotted using the radial velocities of the 
primary and secondary stars, Kw ~ 93 ± 10 kms~^ and 
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Figure 5. Trailed spectra of Hell A4686 A, the nebular lines and Hell A5412 A. The second panel from the left shows the Hell A4686 A 
line computed from the Doppler map. The gaps around phases and 1 in the computed data correspond to eclipse data, which were 
omitted from the fit; the gaps around phases 0.3, 0.5, 1.3 and 1.5 are due to missing data. The data only cover one cycle, but are folded 
over for clarity. 



Kji = 170 ± 10 kms , derived in Section 4.9. The series 



of circles along the gas stream mark the distance from the 
white dwarf at intervals of O.lLi, ranging from l.OI/i at the 
red star to 0.3Li. 

A ring-like emission distribution, characteristic of a Ke- 
plerian accretion disc centred on the white dwarf, is clearly 
seen in Fig. m. There is some evidence for an increase in 
emission downstream from where the computed gas stream 
joins the accretion disc, implying the presence of a bright 
spot. This kind of behaviour has been seen in other CVs - 



e.g. WZ Sge (Spruit & Rutten 1998) 



4.6 Radial velocity of the white dwarf 

The continuum-subtracted data were binned onto a constant 
velocity interval scale about each of the two helium line rest 
wavelengths. In order to m easure the velocities, we use d the 
double- Gaussian method of pchneider fc Young (1980) , since 



this technique is sensitive mainly to the motion of the line 
wings and should therefore reflect the motion of the white 
dwarf with the highest reliability. We varied the Gaussian 
widths between 150-300 kms~^ at 50kms~^ intervals, as 
well as varying their separation a from 200-1500 kms~^. We 
then fitted 



1^ = 7- A'sin[27r(</>-(jf>o)] 



(2) 



to each set of measurements, omitting the 9 points around 
primary eclipse which were affected by the rotational distur- 
bance. An example of a radial velocity curve obtained for 
the Hell A4686A line for a Gaussian width of 300 kms~^ 
and separation 1000 kms~^ is shown in Fig. [tI The data for 
the Hell A5412A line were too noisy to determine accurate 
radial velocities. 

The radial velocity curve has a negligible phase offset, 
an indication that the helium line is a reliable representation 
of the motion of the white dwarf. The results of the radial 
velocity analysis are displayed in the form of a diagnostic 
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Figure 6. Doppler map of Hell A4686A emission in U Sco. 



diagram in Fig. |8|. By plotting K, its associated fractional 
error ax/K, 7 and (po as functions of the Gaussian separa- 
tion, it is possible to select the value of K that most closely 



repr esen ts the actual Kw ( Shaffer, Szkody fc Thorstensen 
198 q). { f the emission were disc dominated, one would ex- 
pect the solution for K to asymptotically reach Kw when 
the Gaussian separation becomes sufficiently large, and fur- 
thermo 



Figure 7. Radial velocity curve of Hell A4686A, measured using 
a double-Gaussian fit with a separation of 1000 kms~^. Points 
marked by open circles were not included in the radial velocity 
fit, as they were affected by the rotational disturbance during pri- 
mary eclipse. The horizontal dashed line represents the systemic 
velocity. 



4.7 Radial velocity of the secondary star 

The secondary star in U Sco is observable through weak ab- 
sorption lines, although many features have been drowned 
out by the nebular lines of the outburst. We compared re- 
gions of the spectra rich in absorption lines with several tem- 
plate dwarfs of spectral types F0-K2, the spectra of which 
are plotted in Fig. llOl The absorption features are too weak 
for the normal technique of cross-correlation to be success- 
ful in finding the value of Kn, but it is possible to exploit 
these features to obtain an estimate of Kr using th e tech- 



e, one would expect 0o 1^> fa.ll LtJ 0. This is seyu Lo 



nique of skew mapping. This t echnique is descibed by smith 



t a separaCiuu of '^1200 km&~ , coiiespondiug Co 
Kw ~ 93 kms~^. There is, however, no sudden increase in 
ax/K, prompting us to employ a difi'erent approach. 

Marsh (1988) suggests that the use of a diagnostic dia- 
gram to evaluate Kw does not account for systematic dis- 
tortion of the radial velocity curve. We therefore attempted 
to make use o f the light centres method, as described by 
Marsh (1988), In the co-rotating co-ordinate system, the 
white dwarf has velocity {0,—Kw), and symmetric emis- 
sion, say from a disc, would be centred at that point. By 
plotting Kx ~ —A' sin 00 versus Ky = —K cos4>o for the dif- 
ferent radial velocity fits (Fig. Bl), one finds that the points 
move closer to the Ky axis with increasing Gaussian sepa- 
ration. A simple distortion which only affects low velocities, 
such as a bright spot, would result in this pattern, equiv- 
alent to a decrease in distortion as one measures emission 
closer to the velocity of the primary star. By extrapolating 
the last point on the light centre diagram to the Ky axis, 
we measure the radial velocity semi-amplitude of the white 
dwarf Kw = 93 ± 10 kms~^. The error on Kw was esti- 
mated from the uncertainty in crossing the Ky axis, given 
the uncertainties associated with the points on the light cen- 
tres diagram. Note the small scale on the a;-axis of the light 
centres diagram. 



Cameron fc Tucknott (1993) 



The first step was to shift the spectral type template 
stars to correct for their radial velocities. We then normal- 
ized each spectrum by dividing by a spline fit to the contin- 
uum and then subtracting 1 to set the continuum to zero. 
This ensures that line strength is preserved along the spec- 
trum. The U Sco spectra were also normalized in the same 
way. The template spectra were artificially broadened by 16 
kms"'^ to account for orbital smearing of the U Sco spectra 
through the 1500 s exposures. The rotational velocity of the 
secondary star was found by the equation 



g(l + g)' = 9.i 



Ki 



(3) 



(^mith fc Dhillon 1998 ) using estimated values of q and Kr 



in the first instance, then iterating to find the best fit values 
given in Section 4.9. The templates were then broadened by 
this value of wsini = 88 kms~^. Regions of the spectrum 
devoid of emission lines (AA5085 — 5395A) were then cross- 
correlated with each of the templates yielding a time series 
of cross-correlation functions (CCFs) for each template star. 
To produce the skew map (Fig. fni), these CCFs were 
back-projected in the same way as time-resolved spectr a in 
standard Doppler tomography ( Marsh fc Home 1988 ). If 



there is a detectable secondary star, we would expect a peak 
at {Q,Kii) in the skew map. This can be repeated for each 
of the templates. 
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Figure 8. The diagnostic diagram for U Sco based on the double- 
Gaussian radial velocity fits to Hell A4686A. 



When we first back-projected the CCFs, the peak in 
each skew map was seen to be displaced to the blue of the 
Kj: = axis by around 30 kms^'^. The reason was that 
we assumed that the centre of mass of the system was at 
rest. The systemic velocity of U Sco derived from the radial 
velocity curve fits in Section 4.6 is 7 = 47 ± 10 kms~^. 
Applying this 7 velocity shifts the peak to the red, towards 
the Kx ^ axis. As Tf^ < 
given by the equation 



Ky, then the value of Kr is 



Kr = (Tff + K 



20/2 



Ky 



SO Kr changes little whilst varying 7. 

The skew maps produced using each of the tem- 
plate stars with 7 = kms~^ show well-defined peaks 
at {Kx,Ky) ~ (—30,160) kms~^. However, applying 7 = 
47kms~^, the peaks in the skew maps shift to {Kx,Ky) ~ 
(15, 180) kms~^. We suspect that the true 7 value of U Sco 
falls between these two limits. The skew map shown in 
Fig. hll is for the K2V spectral type template with 7 = 
47 kms~^, and is arguably the best fitting template for the 
secondary star. The skew map shows a clear peak around 
{Kx,Ky) ~ (10, 165) kms^^. To bring the skew map peak 
to lie on Kx = 0, the 7 value must be modified to ~ 30 
kms~^, which still falls within 2a of the 7 derived from the 
emission lines. The resulting peak falls on Ky = 170 kms~^, 
which we adopt as the radial velocity semi-amplitude for 
the secondary star: Kr — 170 ± 10 kms~^. The uncertainty 
in Kr has been determined from the scatter in Kr due to 
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Figure 9. Light centres diagram for Hell A4686A in U Sco. 

four sources of error: varying the 7 velocity; using different 
spectral type templates; producing skew maps from random 
subsets of the data; the accuracy in measuring the peak in 
the skew map. 

The bottom panel of Fig. y_l| shows a sine wave in the 
trailed CCFs, demonstrating that the peak in the skew map 
is not due to noise in the CCFs. It can be seen that the peaks 
in the CCFs are most prominent around primary eclipse 
(0 ~ 0.7-1.2). There are two possible explanations for this. 
First, by considering the geometry of the secondary eclipse 
in U Sco, we estimate that a maximum of 30 per cent of the 
visible surface of U Sco wil l be obscured by the large ac- 
cretion disc (see Section [l.9| ) between phases 0.35-0.65. We 
would therefore expect the peaks in the CCFs to weaken 
considerably at these phases, which is exactly as observed 
in Figure hll. This explanation is further supported by the 
tentative evidence for a secondary eclipse presented in the 
continuum light curve (Figure W) and it does not affect the 



values for the masses we have derived in Section 4.£ . Second, 



it is possible that the absorption features on the inner hemi- 
sphere of the secondary star are weakened by irradiation due 
to the recent outburst, which would again explain the loss of 
the peaks in the CCFs around phase 0.5. If the latter expla- 
nation were true, it is possible that we have overestimated 
the value of Kr, which would in turn lead us to overestimate 
M\ (see Figure 113). The only way we can reliably account for 
any systematic effects due to irradiation would be to obtain 
higher signal-to-noise observations during quiescence and so 
measure v sin i during primary eclipse. 



4.8 Rotational velocity and spectral type of the 
secondary star 

Mass estimates of CVs using emission line studies suffer from 
systematic errors, and ha ve been thor oughly discussed by a 
number of authors (e.g. 



Marsh 



More reliable mass 
estimates can be obtained by using Kr in conjunction with 
the rotational velocity of the secondary star, usinz. We at- 
tempted to determine v sin i using the procedure outlined in 



Smith, Dhillon & Marsh (1998) , but the data were too noisy 



Previous estimates of t he secondary spectral type have 
ranged from F5V to K2III. [Hanes (1985)| suggests G0±5III- 
V, corroborated by Webbink et al. (1987) with a GUI sub- 
giant. Recent estimates have agreed on the subgiant nature 
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Wavelength (A) 

Figure 10. Normalized spectra of ten template stars. The offset 
between each is 0.035. The U Sco spectrum is a normalized aver- 
age of 43 of the 51 spectra observed; 8 spectra were omitted due 
to cosmic rays interfering with secondary features. The template 
spectra have been broadened by 16 kms~^ to account for orbital 
smearing of the U Sco spectra during exposure and by 88 kms~^ 
to account for the rotational broadening of the lines in U Sco. 



of the secondary. Schaefer (1990)| estimated the spectr al type 



as G3- 6 from the colour s at minimum light, whereas John- 
ston & Kulkarni (1992) measured the observed ratio of cal- 
cium lines to H(5 and found F8±2. Based on the detection 
of the MgI6 absorption in the 1979 outburst spectrum, Ka- 
habka et al. (1999) find the secondary to be consistent with 
a low-mass su bgiant secondary of spectral ty pe K2. This is 
supported by Anupama & Dewangan (2000), who estimate 
the secondary to be a K2 subgiant based on the indices of 
the Mgl& absorption band and the Pel -I- Cal absorption 
feature in the late-decline spectrum of the 1999 outburst. 
Unfortunately, as can be seen in Fig. [lOl it is impossible to 
confirm a secondary spectral type using our data, although 
the MgI6 complex is clearly present. 



4.9 System Parameters 

The radial velocity curves show little or no phase shift, giv- 
ing us confidence that our measurement of Kw reflects the 
motion of the white dwarf. Hence, together with Kr, our 
newly derived period and a measurement of the eclipse full 
width at half-depth (A0i/2), we can proceed to calculate ac- 
curate masses free of many of the assumptions which com- 
monly plague CV mass determinations. Our white dwarf 
radial velocity semi-amplitude, Kw = 93 ± 10 kms~^ com- 
pares favourably with the previous measurement of Kw ~ 
95 ± 50 kms~^ by Barlow et al. (19 81) after the 1979 
outburst. Schaefer fc Ringwald (1995) derived the values 
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Figure 11. The top panel shows the skew map of U Sco for the 
K2V spectral type template. The bottom panel shows the trailed 
CCFs for the same template, clearly demonstrating that the peak 
in the skew map represents Kn and is not an artefact of noise in 
the CCFs. The dashed line shows our adopted parameters for the 
secondary star orbit: Kji = 170 kms~^, 7 = 47 kms"'^ and zero 
phase shift. The data have been folded over for clarity. 



Kw = 87 ± 32 kms"^ and Kr = 5G ± 27 kms"\ but con- 
cluded that these results were unreliable due to inconsistent 
7 velocities, a phase shift in the emission lines and a large 
scatter in the radial velocity curves. 

In order to determine A(j)i^2, we estimated the levels 
of flux outside the eclipse (the principle source of the error) 
and at eclipse minimum, and then measured the full width 
of eclipse half way between these points. The eclipse full 
width at half-depth for U Sco was measured to be A(/)i/2 ~ 
0.104 ± 0.01, in good agreement with A^i/a = 0.106 ± 0.01 
cited by Schaefer fc Ringwald (1995)| . 

W e have opted for a Monte Ca rlo approach similar to 
that of Home, Welsh & Wade (1993) to calculate the system 
parameters and their errors. For a given set of Kw, Kr, 
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Figure 12. Constraints on the masses of the stars in U See. Each dot represents an (Mi, M2) pair. The dashed hnes are hnes of constant 
inclination {i = 80°, i = 85°, i = 90°). The sohd curves satisfy the constraints on Kt^ = 93 ± lOkms"^ and Kji = 170 ± lOkms"^. 



A(;/>i/2 and P, other system parameters are calculated as 
follows. 

The mass ratio, q, can be determined from the ratio of 
the radial velocities 



Ml 
M2 



Kw 

Kr' 



(4) 



7?2/cs can be estimated because we know that the secondary 
star fills its Roche lobe (as there is an accretion disc present 
and hence mass transfer). R2 is the equatorial radius of 
the secondary star and a is the bin ary separation. We used 
Eggleton's formula ( Eggleton 1983 ) which gives the volume- 
equivalent radius of the Roche lobe to better than 1 per cent, 
which is close to the equatorial radius of the secondary star 
as seen during eclipse. 



a 



Q.AQq 



,2/3 



0.6g2/3+ln(l + gi/3)' 



(5) 



By conside ring the geometry of a point ec lipse by a spherical 
body (e.g. Dhillon, Marsh & Jones 1991), the radius of the 
secondary can be shown to be 



-tX.0 \ ^ '? 

— ) = sin 7rA(^i/2 + cos 7rA(;Ai/2 cos i. 



(6) 



Using the value of i?2/ffl obtained from equation H, combined 



with equation |6[ the inclination of the system can be found. 
Kepler's Third Law gives 



K'jPorb _ Ml sin^ 
27VG ' (1 + qf 



(7) 



which, with the values of q and i calculated using equa- 
tions y, y & y, gives the mass of the primary star. The mass 
of the secondary star can then be obtained using equation W. 
The distance of each component to the centre of mass of the 
system (ai,a2), is given by 



0.(1,2) 



K{W,R)P 

2n sin i 



(8) 



The separation of the two components, a, is 01+02, allowing 
the radius of the secondary star, R2, to be calculated from 
equation 0. 

The Monte Carlo simulation takes 10 000 values of Kw, 
Kn, and A(/)i/2 (the error on the period is deemed to be 
negligible in comparison to the errors on Kw, Kr, and 
^01/2), treating each as being normally distributed about 
their measured values with standard deviations equal to the 
errors on the measurements. We then calculate the masses 
of the components, the inclination of the system, the ra- 
dius of the secondary star, and the separation of the compo- 
nents, as outlined above, omitting {Kw,Kii,A<j)i/2) triplets 
which are inconsistent with sini < 1. Each accepted Mi, M2 
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Table 4. System parameters for U Sco 



Parameter 



Measured Value Monte Carlo Value 



Porb (d) 

Kw (kms~^) 
Ku (kms-l) 

9 



1.2305522 
93 ±10 
170 ± 10 
0.104 ±0.010 



95 ±9 

169 ± 10 

0.098 ± 0.007 



0.55 ±0.07 
00.710.0 



ture of RNe, the me chanism depends on the individual sys- 
tem dWebbink 1990|). In the case of U Sco, the TNR model 



was suggested in the light of severe problems with the disc 
instability model (Webbink et al. 1987). The motivation be- 
hind measuring the mass of the white dwarf in U Sco was to 
give positive confirmation that outbursts are the result of a 
TNR, which predicts that the higher the mass of the white 
dwarf, the more frequent the outbursts. For eruptions to re- 
cur on the timescales seen in U Sco (~ 8 years), the mass 
of the white dwarf must be very close to the Chandra s ekhar 
mass of 1.378Af(T, fetarrfield. Sparks fc Truran 1985|; IStar " 



M2/M0 

R2/R0 

a/RQ 

A<f> 

Rd/Ri 



1 55 ±0 7A 



0.88 ±0.17 
2.1 ±0.2 
6.5 ±0.4 



0.12 ±0.01 
0.87 ±0.11 



rficld, Sparks fc Shaviv 1988 Nomoto, Thielemann & Yokoi 
1984|). Our observations allow us to calculate the mass of the 



white dwarf in U Sco to be Afi = 1.55±O.24M0 . Despite the 
large error on the value, we can conclude that U Sco con- 
tains a high mass white dwarf, although additional work, 
such as an accurate v sin i determination during quiescence, 
must be done to further constrain this. 



pair is then plotted as a point in Fig. [13, and the masses 
and their errors are computed from the mean and standard 
deviation of the distribution of spots. The solid curves in 
Fig. h2 satisfy the white dwarf radial velocity constraint, 
Kw = 93 ± lOkms"^ and the secondary star radial ve- 
locity constraint, Kr = 170 ± lOkms"^. We find that 
Ml = 1.55 ± 0.24Mq and M2 = 0.88 ± O.ITM©. The in- 
cMnation of the system is calculated to be i = 82.7° ± 2.9°, 
consistent with the nature of the deep eclipse seen in the 
lightcurves and the presence of double-peaked helium emis- 
sion. 

We computed the radius of the accr etion disc in U Sco 
using the geometric method outlined in Dhillon, Marsh & 



Jones (1991). The phase half-width of eclipse at maximum 
intensity was found to be A(j> = 0.12 ± 0.01 from Figure H 
Combining A^ with q and i derived above produces an ac- 
cretion disc radius of Rd ~ 0.87 ± 0.11i?i, where Ri is the 
volume radius of the primary's Roche lobe. O ur values for 



A(l) an d Rp are co nsistent with those found by Harrop-AUin 
& Wari lier (199"6)1 A^ = 0.129 ± 0.008 an d Rp > 0.92_Ri 



The empirical relation obtained by Smith fc Dhillon 
(1998)|petween mass and radius for the secondary stars in 
CVs is given by. 



R 

Rp> 



(0.93 ±0.09) 




(0.06 ±0.03). 



(9) 



This predicts that if the secondary star is on the main- 
sequence, it should have a radius of O.88i?0. The value we 
find is 2.1 ± 0.2Rq, clearly indicating that the white dwarf 
in U Sco has an evolved companion. 

The values of all the system parameters of U Sco derived 
in this section are listed in Table 0. 



5 DISCUSSION 

5.1 TNR Model 

The outburst events in RNe could be the result of the TNR 
model, as seen in classical novae, or the disc instability model 
of dwarf novae. It is likely that given the heterogeneous na- 



5.2 U Sco as a Type la SN Progenitor 

Type la supernovae (SNe la) are widely believed to be ther- 
monuclear explosions of mass-accretin g white dwarfs (see 
Nomoto, Iwamoto & Kishimoto 1997 for a review). The 



favoured model for progenitor binary systems is the Chan- 
drasekhar mass model, in which a mass-accreting carbon- 
oxygen white dwarf grows in mass up to the Chandrasekar 
limit and then explodes as a SN la. For the evolution of 
accreting WDs towards the Chandrasekhar mass, two sce- 
narios have been proposed. The first is a double degener- 
ate scenario, where two white dwarfs merge to cross the 
Chandrasekhar limit. A candidate for this has been iden- 



tified in KPD 1930+2752 ( [Maxted, Marsh fc North 2000| ). 
The second is a single degenerate scenario, where matter 
is accreted via mass transfer from a stellar companion. A 
new evolutionary path for single degenerate progenitor sys- 
tems is discussed by Hachisu et al. (1999), which descibes 
how the secondary (a slightly evolved main-sequence star) 
becomes helium rich. Their progenitor model predicts he- 
lium enriched matter accretion onto a white dwarf. Objects 
which display these characteristics in the form of strong 
Hell A4686A lines are luminous super-soft X-ray sources 
and certain recurrent novae, such as U Sco and the others in 
the U Sco subclass (V394 CrA and LMC-RN). The obser- 
vations in this paper, in confirming the helium rich accreted 
matter, and constraining the white dwarf mass to be near 
the Chandrasekhar mass, make U Sco the best candidate for 
this evolutionary path to a SN la. 

The time to supernova can be estimated by dividing 
the amount of mass required to reach the Chandrasekhar 
mass by the average mass accretion rate (Mav)- Our results 
suggest a minimum white dwarf mass of --^ I.3IM0, which 
requires it to accrete approximately O.O7M0 to become a 
SN la. In their model, Hachisu et al. (2000b) find that the 
white dwarf in U Sco grows in mass at an average rate of 
Mav ~ 1.0 X lO~^Af0 y~^. Using this value, we therefore 
predict U Sco to become a supernova within ~ 700 000 years. 



6 CONCLUSIONS 

(i) We have shown that U Sco contains a high mass white 
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dwarf (Ml = 1.55 ± O.24M0), confirming tliat tlie outburst 
mechanism is the TNR model. 

(ii) The presence of an accretion disc has been confirmed 
for the first time by clear evidence of rotational disturbance 
in the Hell emission lines during eclipse, as well as their 
double-peaked nature. 

(iii) The secondary star has a radius of 2.1 ± 0.2Rq, con- 
sistent with the idea that it is evolved. 

(iv) The mass of the white dwarf in U Sco is Mi — 1.55 ± 
O.24M0, implying that it is the best SN la candidate known 
and is expected to explode in ~ 700 000 years. 
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